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ABSTRACT 

Very high-energy 7-rays (VHE; E>100 GeV) have been detected from the direction of the 
Galactic Centre up to energies E>10 TeV Up to now, the origin of this emission is unknown 
due to the limited positional accuracy of the observing instruments. One of the counterpart 
candidates is the super-massive black hole (SMBH) Sgr A*. If the VHE emission is produced 
within w 10 15 cm « 1000 r G (r G = GM/c 2 is the Schwarzschild radius) of the SMBH, 
a decrease of the VHE photon flux in the energy range 100-300 GeV is expected whenever 
an early type or giant star approaches the line of sight within « milli-arcseconds (mas). The 
dimming of the flux is due to absorption by pair-production of the VHE photons in the soft 
photon field of the star, an effect we refer to as pair-production eclipse (PPE). Based upon the 
currently known orbits of stars in the inner arcsecond of the Galaxy we find that PPEs lead 
to a systematic dimming in the 100-300 GeV band at the level of a few per cent and lasts for 
several weeks. Since the PPE affects only a narrow energy band and is well correlated with 
the passage of the star, it can be clearly discriminated against other systematic or even source- 
intrinsic effects. While the effect is too small to be observable with the current generation 
of VHE detectors, upcoming high count-rate experiments like the Cherenkov telescope array 
(CTA) will be sufficiently sensitive. Measuring the temporal signature of the PPE bears the 
potential to locate the position and size of the VHE emitting region within the inner 1000 tq 
or in the case of a non-detection exclude the immediate environment of the SMBH as the site 
of 7-ray production altogether. 
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1 INTRODUCTION 



The central 



region of our Galaxy harbours a 
super-massive black hole (SMBH) with a mass 
M = (4.31 ± 0.06| s tat ±0.36| Ro ) x 10 6 M o at a distanc e 
from Earth of R = 8.33 ± 0.35 kpc (IGillessen etaD l2009h . 
Its proximity enables detailed observations with high angular 
resolution in the radio and, recently by the use of adaptive optics , 
in the near infrared band (for a review see e.g. iReidl [2009). 
Accurate measurement of the Kepl erian orbits of stars in the inner 
arcse c ond of the Galactic Ce n tre dSchodel et al.ll2002l ; iGhez et al.l 
120031 : lEisenhauer etaD 120051 : iGillessen et aD l2009h has revealed 
the presence of a co mpact massive object co-located with th e 
radio-source Sgr A* tealick & Brownl 1 1974 iReid et all [200% 
The astrometry allows for a shift of up to 2 mas (IGillessen et al.l 
2009) being the upper limit on the systematic error of the position 
between the radio and near infrared positions of Sgr A* . 

Recently, very high-energy (VHE; E>100 Ge V) emission has 



been d etected from t he Ga la ctic Centre reg ion (Aharonian et al 
2004; iKosack et al] |2004 iTsuchiva et ail |2004 lAlbert et al 



2006). A number of different scenarios for particle accelera- 



tion and 7-ray produc t ion ha ve been discussed in the literature. 
Ah aronian &Neronov (2005) have considered 7-ray production 
both from electrons via curvature radiation as well as from hadronic 
processes (photo-meson production and proton-proton scattering). 
The energetic photons can be produced as close as 20 tg with- 
out suffering strong absorption via pair-production with the photon 
field generated by the low-lu minosity SMBH (L < 10~ 8 L E dd) 
(I Aharonian & Neronovll2005|). In anothe r SMBH-related scenario 
suggested by lAtovan & D ermer (2004) an assumed relativistic 
wind injected in the vicinity of the SMBH would lead to the for- 
mation of a synchrotron/inverse-Compton nebula similar to pulsar 
wind generated nebulae. In this case the VHE emission would be 
produced in the vicinity of the shock that is leading to accelera- 
tion of particles at a distance of (D(10 16 ) cm to the SMBH. More 
production scenarios, which are not directly related to the SMBH, 
have been discussed: T he pulsar-wind nebula candi date G359.95- 
0.04 dWang e7aT]l2006h . the stellar cluster IRS 13 iMaillard et all 
120041) . the low-mass X-ray binary candidate J174540.0-290031, 
Sgr A East ( Aharonian et al. 2006a) or self- annihilating Dark Mat- 
ter (see e.g. lHornsl2005g) . Fig. Q] shows an X-ray image (0.3-8 keV) 
derived from 156 hours of archival Chandra observations including 
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Figure 2. A cartoon of the line of sight integration and the used geometry. 




Figure 1. Chandra X-ray image of the Galactic Centre region accu- 
mulated from 156 hours of data in the 0.3-8 keV band. The possible 
sources for the VHE emission are marked, the small circle indicates 
the position of the H.E.S.S. source (I = 359 °56^4L1^ ± 6.4" (stat.), 
b = -0°2 / 39.2 // ± 5.9" (stat.)) as reported in lAcero et alJ (2009). The 
large white ellipse marks the combined statistical and systematical (6 arc- 
sec) H.E.S.S. positional uncertainty region. The supernova candidate Sgr A 
East is outside of the image. 



the X-ray emission observed from the various counterpart candi- 
dates located within the location uncertainty of the VHE source. 

Re cent improved rec onstruction of the location of the VHE- 
source (lAcero et al.l 2009) has excluded the emission originating 
predominantly from Sgr A East, while the other possibilities men- 
tioned above remain viable. 

The absence of flux variability in the VHE band during simul- 
taneous observations of Sgr A* by H.E.S.S. (High Energy Stereo- 
scopic Systme) and Chandra seem to disfavour a common origin 
of the X-r ay and the VHE emis sion, yet such a scenario can not be 
excluded ( Aharonian et al.ll2008h . 

Despite the efforts to reduce the positional error of the VHE 
source, the current uncertainty of 6 arcsec (syst.) per telescope axis 
(lAcero et al.ll2009h is limited to systematic errors. The prospects 
for a substantial improvement on the systematic uncertainties are 
limited by the mechanical stability of air Cherenkov telescopes, 
which would require substantial improvement with respect to the 
currently used light-weight construction. The Galactic Centre re- 
gion is too crowded to identify the VHE source reliably. Unless a 
clear timing signature can be established, this will not change. 

In this paper we propose the use of a time-dependent attenu- 
ation feature, which provides a novel opportunity to constrain the 
size and location of the Galactic Centre VHE- source with an ac- 
curacy of milli- arc seconds in the case of the VHE emission being 
produced in the vicinity of Sgr A* . We provide a long-term (60 yrs) 
prediction of flux variations in different energy bands which can be 
caused by the passage of stars close to the line of sight towards the 
VHE source. 



2 PAIR-PRODUCTION ECLIPSES 

In the following we assume that VHE 7-rays originate from a re- 
gion near the SMBH Sgr A* . These VHE photons interact with 
soft photon fields via the pair-production process (7vhe + 7* — >• 
e + + e~). Stars on stable Keplerian orbits in the direct vicinity of 



Sgr A* (< 1") provide these (low-energy) photon fields, and pair- 
production can lead to an energy-dependent dimming of the VHE- 
flux. Since the position of the star changes, the optical depth r for 
the VHE photon varies along the orbit of the star. In the following 
we will calculate r as a function of time and 7-ray energy for stars 
on Keplerian orbits. 

2.1 Pair-production 

The total cross-section for pair-production (after averaging over po- 
larisation states) is given as 

3<7t , 
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with gt being the Thomson cross-section, f3 — (1 — 1/s) 1 ^ 2 
and s — £E/(2m 2 e c 4 )(l — cos#) the (normalized) centre-of- 
momentum energy of two photons of energy e and E interac ting at 
an angle 6 in the laboratory frame (iGould & SchrederJ [l967)Fl We 
choose the laboratory frame as the system where the VHE source 
is at rest. The cross-section is strongly energy-dependent and has a 
maximum for photon energies of 

£ max = .9(^) -1 TeV, (2) 

corresponding to a maximal cross-section for a 1 TeV 7-photon 
interacting with a 0.9 eV (A w 1.38 /im, near infrared) stellar pho- 
ton. The relative velocities of the observer as well as the star with 
respect to the VHE source are considered to be small (with respect 
to the speed of light c) and will be neglected in the following. For 
the considered case we assume the stellar atmosphere of a star with 
radius i?* and effective temperature T e ff to be the source of low- 
energy photons. The differential photon density dn = n(e) de at 
the photosphere of the star can be approximated with a black-body 
spectrum: 

h 3 c 3 



n(e) 



exp 



(3) 



with h being Planck's constant, c the speed of light and ks Boltz- 
mann's constant. The optical depth for VHE photons with energy 
E is calculated by integrating along the line of sight subtending a 
perpendicular distance D± to the star (see Fig.|2]for a definition of 
the geometry): 

Rq 1 2"7T OO 

t(E) = J dl J dcosG J dcj) J dea(s,E,0)n(e) (l-cos0), 

c min £ min 

(4) 

1 In this paper e denotes the energy of the low-energy stellar photon, E 
is the energy of the high-energy photon coming from Sgr A* ; m e is the 
electron mass. 
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Figure 3. Energy dependence of the absorption exp (— r) for a fixed per- 
pendicular angular separation of 1 mas (for a distance Rq = 8.33 kpc). The 
three curves show the absorption for temperatures T e ff of the star; the stellar 
radius is adjusted according to the range of stellar types suggested for the 
well-studied S2 star (T eff = 19 000 K, R* = 13.1 R®;T^ = 25 000 K, 
= 10.7 R®; T eff = 30 000 K, R* = 9.4 R®) CSartins et all2008h . 



with c min = [l - (tt ) 2 ] V2 and e min = 2m 2 e c 4 /[E(l - cosfl)] 
being the minimum energy for pair-production to occur (e.g jDubusl 
2006). For R* <C D the integration over the solid angle in Eq.01 



leads to the point- source approximation: 



R 



(E) = J dl J den^j^j a(e,E,0)n(e) (l-cos0). (5) 



For the distances examined here this is an appropriate approxima- 
tion and will be used in the following. In order to demonstrate the 
uncertainty on the optical depth that results from the uncertainty 
on the stellar parameters, we calculate the optical depth for the 
well-studied star "S2". This star is constrained to have a temper- 
ature 19 000 K < Teff < 3 000 K and theref ore stellar radii 
13.1 R® > R* > 9.4 R® (Mar tins etall feoOS) corresponding 
to spectral types B0-B2.5V. Fig. [3] shows the result of the calcu- 
lation of Eq. \5\ for a fixed distance D± corresponding to a pro- 
jected angular separation of 1 mas of the VHE source and the star. 
The absorption shows a pronounced minimum around 200 GeV for 
T e ff = 19 000 K shifting to smaller energies (see also Eq.[3 with 
increasing stellar temperature. The stellar types have been deter- 
mined spectroscopically for nearly 50 so-called S-stars. Most stars 
are early type, B -main- sequence- stars with s tellar radii of 10-2 
R® and effective temperatures of 2-3 x 10 4 K ( Mar tins et al.l2 008). 

As shown in Appendix [A| an additional contribution of cas- 
cade radiation coming from the electron/positron pairs produced in 
the PPE is negligible. The resulting change is similar to the un- 
certainties resulting from the uncertainties of the orbital elements. 
Therefore, the effect is generally not of importance. 

Another effect which could change the timing signature is a 
possible change of the accretion flow towards the SMBH induced 
by the star closely passing Sgr A* . Direct wind-driven accretion 
can lead to an increased accretion rate, but the fast motion of the S- 
stars does not allow any accretion during periastron passage jLoebl 
120041) . Furthermore, the PPE can be treated as independent from ac- 
cretion induced variability, because it is strongly time- and energy- 
dependent. 
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Figure 4. The projected (on the plane of the sky) orbits of three selected 
S-stars (see text). The cross marks the position of Sgr A* . 
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Figure 5. Absorption exp[— r(E, £)] due to pair-production eclipses for 
three S-stars for VHE photons of energy E = 200 GeV as function of the 
calendar year. 



2.2 Stellar orbits and optical depth for VHE photons 

The temporal evolution of the optical depth calculated by means of 
Eq.[5]depends on the position of the star relative to the line of sight, 
which in turn varies along the orbit of the star. The orbital parame- 
ters of 29 stars in the central arcsecond of the Milky Way have been 
well-measured using both spectroscopy and time-dependent as- 
trometry with unprecedented accuracy utilising the VLT advanced 
adaptive optics in the ne ar-infrared (^T-)band feisenhauer et al.l 
l2005l : lGillessen et alJ2009h . 

With the given orbital parameters we derived the three- 
dimensional position usin g the known methods described e.g. in 
Murr ay & Dermottl dl999h . The projection of three stellar orbits on 
the plane of the sky is shown in Fig. |4] The stars S2, S6, and S14 
have been selected in particular, because they approach the line of 
sight closest and therefore are expected to lead to the largest ab- 
sorption effect. Since the spectral types of these stars are not ex- 
actly known, representative values for a B-type main sequence star 
(similar to the parameters for S2) are used: R* = 10.7 R® and 
T eff = 2.5 x 10 4 K. 

Table Q]lists the known orbital parameters, their uncertainties, 
the derived minimal angular separation $ « D±/Ro as well as the 
correspondingly calculated absorption exp(— r) for VHE photons 
of energy E = 200 GeV. The optical depth takes into account the 
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Table 1. Parameters of selected S-stars: The first 7 columns show the orbital parameters taken from Gillessen et al. ( 2009), the last two columns list the minimal 
angular separation of the star to the line of sight # (in mas) and the resulting maximum absorption exp [— r(E)] for E = 200 GeV. 



Star 
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n[°] 
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t P [yr] 
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exp ( — t) 


S2 
S6 
S14 


0.123 ± 0.001 
0.436 ± 0.153 
0.256 ± 0.010 


0.880 ± 0.003 
0.886 ± 0.026 
0.963 ± 0.006 


135.25 ± 0.47 
86.44 ± 0.59 
99.4 ± 1.0 


225.39 ± 0.84 
83.46 ± 0.69 
227.74 ± 0.70 
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Figure 6. Absorption exp[— r(E, t)] due to pair-production eclipses for the 
star S14: the curves denote the most likely absorption and the 99 per cent 
c.l. region. The underlying histogramme indicates the smoothed probability 
density function resulting from the Monte-Carlo calculations (see text for 
details). 



full three-dimensional geometry of the problem. For the three stars 
considered here the calculated maximum attenuation of the VHE- 
flux turns out to be 0.5 per cent for S2, 2.4 per cent for S14 and 
3.4 per cent for S6. Fig. \5\ depicts the predicted absorption of these 
three stars until the year 2070. 

In order to propagate the uncertainties on the orbital parame- 
ters into uncertainties on the absorption, a dedicated Monte-Carlo- 
type simulation has been performed. Since the errors on the dif- 
ferent orbital parameters are correlated, we transformed the pa- 
rameters in the system of the eigenvectors of the covariance ma- 
trix and varied these transformed (then uncorrelated) parameters 
randomly 10000 times following a normal distribution within the 
known variance. After back- transforming, we end up with 10000 
orbital parameter- sets tracing the probability density function in- 
cluding the correlation of the parameters. 

In order to derive the resulting uncertainty on the absorp- 
tion effect, we calculated exp[— r(E, £)] for each parameter-set ob- 
tained in the Monte-Carlo simulation. As a consequence of the un- 
certainties of the orbital parameters the time when the minimum 
of exp(— r) occurs varies considerably. Since we are mainly in- 
terested in the uncertainty on the absorption at a fixed time we 
shifted all curves in time to the mean (best-fitting) value. The most 
likely attenuation exp(— r) for the star S14 is shown in Fig. [^to- 
gether with the 99 per cent confidence interval. For the other stars 
the uncertainties (68 per cent confidence level) on the minimum of 
exp(— t) are listed in Table [T] Moreover, the often unknown ex- 
act stellar type leads to additional uncertainties on the absorption, 
which are of similar magnitude as the uncertainties due to the errors 
on the orbital elements. 
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Figure 7. Effective areas for different IACTs used in the simulation of the 
light curves. 



2.3 Simulated light curves 

The predicted flux modulation due to PPE is at the level of a few per 
cent (depending on the estimated uncertainties of the orbit and stel- 
lar type). In order to investigate the possibility of discovering PPEs 
measurements of light curves are simulated for future Imaging Air 
Cherenkov Telescopes (IACTs). For completene ss and compa rison 
light curves for the H.E.S.S. telescope system fantonll2004l) are 
simulated as well. 

The expected photon rate R(Ei,E2,t) in an energy band 
E\-E<i is calculated: 

E 2 

R(E 1 ,E 2 ,t) = J dE A eS (E)<j>(E) exp[-r(E,t)] (6) 

with A e ft(E) the effective area of the IACT, <j>{E) the differen- 
tial VHE photon spectrum from the Galactic Centre, and r(E,t) 
the optical depth calculated by means of Eq. \5\ The photon spec- 
trum <j>( E) is assumed to follow a power-law with photon in- 
dex 2.25 dAharonian et al.ll2006ah . The recently published finding 
of an exponential cut- off in the energy spectrum at « 16 TeV 
dAharonian et al J [2009b changes the predicted rate between 1 and 
10 TeV only marginally and is therefore neglected here. 

The effec tive areas assumed (see Fig. are parameterised 
based upon lAharonian et alj d2006bh for H.E.S.S. and IPunchl 
d2005h for H.E.S.S. II (see Appendix [0 for further details). For a 
next generation teles cope- system (like the planned CTA, see e.g. 
iHermann et alJ2008h we generically assumed an increase of the ef- 
fective area by a factor of 25 with respect to the H.E.S.S. collection 
area as well as a lower energy threshold of about 10 GeV. In the 
following three separate energy bands are considered: the very low 
energy (10-50 GeV) band, the low energy band (100-300 GeV), 
and the high energy band (1-10 TeV). Given the energy dependence 
of the PPE, most of the attenuation will affect the low energy band, 
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Figure 8. Simulated light-curves for the star S 14 for 3 different experiments. The red curve shows the calculated absorption for E = 200 GeV, the points 
show the normalised photon-rate in the low energy band (100-300 GeV) where the effect has its maximum. The shaded area marks the photon-rate and its 
statistical uncertainty without absorption in the combined high (1-10 TeV) and very low (10-50 GeV) energy band. In this band the expected absorption is 
negligible (see Fig. 0. From left to right: H.E.S.S., H.E.S.S. II, CTA (note: change in scale of the three panels y-axis). 



whereas the adjacent energy bands remain without attenuation. The 
available observation time is calculated for a site at the Tropic of 
Capricorn taking into account the dark (moon-less) visibility time 
of the Galactic Centre between roughly March to September. For 
each new moon period, the number of detected photons is estimated 
in the three energy bands by calculating N p h(Ei, E2) — R t v - ls , 
with tvis being the actual observation time available until the Galac- 
tic Centre reaches an altitude smaller than 30°. The variance on 
the number of observed photons is estimated under the assump- 
tion of Poissonian statistics, and the expected value iV p h is sam- 
pled following a normal distribution. This is done separately for 
the three energy bands. The resulting simulated light curve for the 
passage of S14 in the year 2047 (normalised to the unattenuated 
light curve) for the low-energy band for the three IACT setups is 
shown in the three panels of Fig. [8] In addition to the attenuated 
low-energy band, the other two energy bands are combined into a 
comparison measurement marked by the shaded region. This com- 
parison measurement allows the elimination of systematic effects 
due to observational artefacts or even source intrinsic, broad band 
variability of the source. It should be noted that the variation of the 
collection area for different zenith angles has not been included in 
the calculation. In principle, for smaller altitude angles, the lowest 
energy band will not be accessible any more. However, this will be 
compensated by an increase of the collection area at larger energies 
not altering the conclusions drawn here. 

The large uncertainty of the time of the closest approach to the 
line of sight for S14 3 yrs) leaves us freedom to shift the light- 
curve in time. The time of the closest approach is shifted slightly 
with respect to the best-fit value by 2 months to the middle of the 
year 2047. In this way the maximum observation time available 
for Cherenkov telescopes is in phase with the time of the closest 
approach. Additionally, we assumed the largest possible attenua- 
tion of the flux within the orbital uncertainties (see Table Q]). Fig. 
[8] shows that a state of the art IACT like H.E.S.S. (left graph) is 
not able to detect a PPE within the measurement errors. Under 
the favourable assumptions used H.E.S.S. II (middle graph) will be 
able to discover a PPE within the errors. A next generation IACT- 
system like the planned CTA (right graph) will discover PPEs even 
under unfavourable conditions. 



ming of the 7-ray flux from the direction of the Galactic Centre. 
Here, we have calculated for the first time the effect of PPEs for the 
Galactic Centre and provide predictions for the next few decades. 

For the currently known orbits of S- stars the effect will be 
detectable with future IACTs. The next close approach to the line 
of sight causing a dimming > 1 per cent will be around the year 
2047. At the moment, S- stars to a limiting magnitude of = 
18 (K-bmd) are trac ked within 2" around Sgr A* with the VLT 
dGilles sen et alj|2009k resulting in a total number of 109 stars rou- 
tinely followed. It is therefore expected that the orbital parameters 
of these, and even fainter stars, will be determined in the future; 
some of these stars may lead to even stronger PPEs at an earlier 
time. 

So far, we have neither considered the effect of binary systems 
nor of non main sequence (e.g. giant) stars leading to a possibly 
larger absorption effect. 

A detection of a PPE has the potential to resolve the current 
source confusion at the Galactic Centre. It could constrain the 7-ray 
emitting region to the direct vicinity of the SMBH Sgr A* within 
roughly 1000 re excluding the other possible VHE emitters in 
the region (the pulsar- wind nebula candidate, the stellar cluster or 
the low-mass X-ray binary system). Moreover, the position of a 
VHE emitting region could be reconstructed with milli- arc second 
accuracy for the first time. A non-detection could exclude the inner 
~ 1000 tg around Sgr A* to be the dominant source of the 7- 
emission from the direction of the Galactic Centre. 

The effect of PPEs can therefore be used to probe the location 
of particle acceleration processes in the environment of a quiescent 
super-massive black hole. 
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3 SUMMARY AND DISCUSSION 

Under the assumption that a compact emission region (< 1000 tg) 
in the vicinity of Sgr A* is the dominant source of the 7-rays, pair- 
production eclipses can lead to a time- and energy-dependent dim- 
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APPENDIX A: RE-SCATTERED PHOTONS FROM 
INVERSE-COMPTON CASCADES 

In the following we calculate the re-scattered flux /J asc , which is 
observed in addition to the attenuated 7-ray flux from Sgr A* : 



nobs 
Jv 



^exp(-r) + 



(Al) 



where j v labels the initial flux from Sgr A* . Since we are mainly 
concerned with the effect of the re- scattered photons on the pair- 
production eclipse, the calculation will only consider the first gen- 
eration of cascade photons produced after pair-production and in- 
verse Compton up- scattering. While the cascade may continue, the 
subsequent scatterings will shift the mean energy of emerging pho- 
tons well below the energies where the pair-eclipse can be observed 
until eventually the energy threshold for pair production is not sur- 
passed any more. The remaining electrons and positrons will con- 
tinue to cool. For reasons of simplified notation, we use in the 
following energy in units of rest mass of electron. Adapting the 
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Figure Al. Geometry used for the calculation of the observed flux of re- 
scattered photons. 



notation in the main body of the text for the high energy photon 
E — emieC 2 , for the soft (stellar) photon e — €2m e c 2 , and for the 
electron/positrons E e = ^m e c 2 . 

Given the fact, that the mean free path length for inverse 
Compton scattering (G(> 10 16 ) cm) is considerably larger than 
the gyro-radius in /xG magnetic fields (G(< 10 14 ) cm), the elec- 
tron/positron pair is quickly isotropised, so that the cascade acts 
like an optically thin diffusor and wave-length shifter, replacing the 
incoming photon of energy ei by two isotropically emitted pho- 
tons of lower energy. One of the two re-emitted photons produced 
by inverse Compton scattering of the leading electron/positron car- 
ries a substantial fraction of the energy of the incoming photon. 
We choose to model the radiation transport by assuming that the 
7-rays are emitted isotropically from the 7-ray photosphere around 
the SMBH, which we define to be the region where r ~ 1 corre- 
sponding to the surface of a sphere of radius « 20 tg- During the 
eclipse a fraction 1 — exp[— r(ei)] of photons of energy ei is ab- 
sorbed along the line of sight subtending an angle a with respect to 
the line connecting the SMBH and the star. The pairs produced are 
then isotropised before upscattering photons from the star to lower 
energies e[ . The effect of cascading is thus a shift of energy from 
ei to e[. Quantitatively, we calculate the probability P(e[ \e±) for a 
given photon of energy a to be re-emitted as a photon of energy e[ . 
The geometry of the system can be simplified by assuming that the 
cascade takes place at the point where the line of sight is closest to 
the stellar surface. This simplification is justified as the probability 
for pair production is highest where the photon density is largest. 
In order to estimate the fraction of the re- scattered flux we also as- 
sume the SMBH as a point source. The resulting re-scattered flux 
can therefore be written in the following way (see Fig. IA1I for a 
sketch of the underlying geometry): 



/r sc (ei) = J de 1 P(e 1 \e 1 )f u (e 1 



da R 2 sin a — 



• exp [-r(a\ a) 



(A2) 
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where the integration over the sphere R = D*/2 incorporates 
all isotropically re-scattered components of the emission from the 
SMBH as seen by the observer. To account for the energy shift 
the integral has to be convolved with P(ei |ei). Integrating over 
and using the geometric relations D' = D* cos a 7 and a = a/2 



Locating the VHE source in the Galactic Centre 7 



reveals 



/rVi)=i J de 1 P(e' 1 \e 1 )f v (e 1 ) 

a max 

x J dot tan a (l — exp [— r(a, ei)] ) 



(A3) 



The limit a m in is given by the radius of the star while awx is 
defined by the 7-photosphere of the SMBH. The geometric effect 
results in a factor of « 5 x 10" 3 for VHE photons of energy 
200 GeV. 

In order to evaluate the re-distribution kernel P(e[ \e± ), we use 
the differential pair-production cross section da/dj in the kine- 
matically accessible range: ei(l — y/l — (e^) -1 ) < 27 < 
e\ (1+ \/l — (ei 62 ) ~ 1 ) - This c ross section is conveniently approx- 
imated (lAg; 
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This cross-section is accurate to the level of a few per cent in the 
limiting case of e\ ^> €2 (see also discussion by Boettcher et al. 
419971)). The assumption of isotropy underlying this approximation 
is not entirely justified. However, when evaluating the total cross 
section and the corresponding optical depth for the chosen geom- 
etry, the relative differences to the line of sight integration of the 
angle dependent cross section is a few per cent only. In order to 
evaluate P, the differential cross section is normalised such that 
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taking into account that a pair of particles is produced. The 
inverse Compton scattering cross section has been taken from 
Blumenth aT& Gouldl dl97oh 
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with z = ei/7 and b = 4e 2 7. Again, this differential cross section 
is normalised, 



dn 



de[ 



7( e i> e 2,7) 



da 



7( e i' e 2,7) 

Jde[^(e' u e, 7 ) 



de[ 
J den(e) 



(Al) 



so that finally 

P(ei|ei) = J dj J den(e)^(ei,e,7) J de n(e)^(e u e\ 7). 

(A8) 

In Fig. lA2l we display the effect of the cascade by showing the 
ratio of the observed flux to the initial flux with and without the ad- 
ditional cascade component. The result indicates a small increase of 
the flux when including the first generation cascade. The resulting 
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Figure A2. Flux ratio fS hs /fu for the star S14 at the time of the largest 
absorption effect with (dashed blue) and without (solid red) the contribution 
of the photons from the cascade. 



Table Bl. Parameters used for the calculation of the effective area for the 
three different IACTs. 





H.E.S.S. 


H.E.S.S. II 


CTA 


91 


6.85 x 10 9 cm 2 


2.05 x 10 10 cm 2 


1.71 x 10 11 cm 2 


92 


0.0891 


0.0891 


0.0891 


93 


5 x 10 5 MeV 


1 x 10 5 MeV 


1 x 10 5 MeV 



change is similar to the uncertainties resulting from the uncertain- 
ties of the orbital elements. The slight bump due to the isotropic 
cascade radiation could in principle be observed even during the 
passage of a star far away from the line of sight (thus leading to no 
absorption effect). The effect is however at the level of 0.1 per cent 
and therefore not easily observable. 

The additional feature (pile-up) in the energy spectrum at ener- 
gies below the main absorption peak has a timing signature shifted 
with respect to the time-dependent absorption due to the cooling 
time of the electrons, which is in O(20) d. The resulting (time- 
dependent) emission spectrum of electrons cooling in the Klein- 
Nishina regime could in principle be c alculated using the approach 
suggested by Manolakou et all d20Q7h . Given the smallness of the 
effect however, we consider it of little interest here. 



APPENDIX B: PARAMETRISATION FOR THE 
EFFECTIVE AREAS 

For the effective areas A e ^(E) used in the simulation of the light- 
curves (Eq. [6) we used the following parametrisation: 



A efi (E) =01 



— V 

MeV J 



exp 



(-§)• 



(Bl) 



The parameters for the three different IACTs are given in Table lBll 
For H.E.S. S. and H.E.S.S. II the c urve is based upon the effective 
area from lAharonian et all j2006bh and lPunchl d2005h respectively. 
For CTA we assume the effective area to increase by a factor of 25 
as given m Hermann et al. (2(X)|). 



